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Abstract. Let X^ngN'' f^^" be a multivariate generating function that converges in a 
neighborhood of the origin of C*. We present a new, muhivariate method for computing 
the asymptotics of the diagonal coefficients fain,...,ajn and show its superiority over the 
standard, univariate diagonal method. 



1. Introduction 

Let -F(x) = X^ngp^d /nx" be a complex power series that converges in a neighborhood 
of the origin but not on all of C^. Here x = {xi, . . . , x^), n = (ni, . . . , na), and /nx" = 

We wish to compute asymptotics for the diagonal coefficients fain,...,adn for fixed positive 
integers ai,...,ad, a task often useful in enumerative combinatorics. For simplicity of 
presentation, we suppose that /n > and that F is rational, although much greater 
generality is possible. Thus far the only general method to extract diagonal asymptotics 
to be found in the literature is what we call the standard diagonal method. It consists of 
two steps: (1) find a closed form or defining equation for the diagonal generating function 
G{x) := ^ngpj fain,...,aanX^ via contour integration (in the style of [Fur 6 7] and [HK71] and 
summarized in |Sta99l Section 6.3]); (2) apply univariate singularity analysis to G{x) (in 
the style of |FSj for instance) to compute the asymptotics. However, as elaborated below, 
this method is quite limited, working well only for main diagonals in two variables, that 
is, for computing the asymptotics of 

In this article we present a new, more versatile diagonal method. It consists of one step: 
apply multivariate singularity analysis (in the style of |PW02| . |P W04] . [BP] ) directly to 
-P(x) to compute the asymptotics. 

2. Limitations of the Standard Diagonal Method 

What is wrong with the standard diagonal method? It certainly works well for computing 
main diagonal asymptotics in two variables, as witnessed by the following example taken 
from |MZSn4j . 

Example 2.1 {Zigzag-free Binary Words). The bivariate generating function 

m,n a > ti a 

counts the number of words over a binary alphabet, {0, 1} say, that have m zeros and n 
ones and do not contain zigzags, that is, the subwords 010 and 101. The main diagonal 
coefficients /„„, then, count zigzag- free binary words with an equal number of zeros and 
ones. 

To compute the asymptotics of /„,„ using the standard diagonal method we proceed 
as follows. Since F{x, y) is rational and holomorphic in a neighborhood of the origin, for 
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fixed X small enough F{x/t, t) will be rational and holomorphic as a function of t in some 
annulus about t = 0. Thus in that annulus it can be represented by a Laurent series whose 
constant term is [t^]F{x/t, t) = J2n>o fnnX^, the series we want. By Cauchy's Integral and 
Residue Theorems we have that for some circle 'jx about t = 



J2fnnX"'=[t^]F{x/t,t) 



2m 4 t 

^Res{F{x/t,t)/t;t = Sk) 

k 



where the are the "small" singularities of F{x/t,t)/t, that is, the ones satisfying 
liuix^o Sk{x) = 0. Since F is rational, these singularities are poles and algebraic func- 
tions of X, so that the residue sum, the diagonal generating function, is also an algebraic 
function of x. 
In particular, 

, I \ I 1 I X I X 

Fix t,t) t = — ^ ^- , 

which has a single simple pole approaching zero as x approaches 0, namely, s = + x — 
x"^ — y/\ — 2x — x"^ — 2x^ + x^). The residue, the diagonal generating function, is then 



, . , x^ +X + 1 

G{x) :-- 



^ - 3x + 1 

Now, the singularity of G closest to the origin is c<j := (3— \/5) /2, and its reciprocal is the 
exponential growth order of the coefficients of G. To determine the leading subexponential 
factor we note that 



/ 



G{x) ~ (1 - -) 

UJ 



^ ^-1/2 



-UJ X 




+ X + 1 1 ^-,-1/2 



V5 



as X — > (jj, so that 

2 nV2-i 



fnn = ffn ~ 



y7fr(i/2) Vs-Vsy 77!;^ 

by asymptotic transfer (see [FSt Chapter VI] for instance). < 

However, the standard diagonal method encounters major problems off the main diago- 
nal even in two variables, as illustrated in the next example adapted from |Sta99t Section 
6.3]. 

Example 2.2 [Lattice Paths). The bivariate generating function 

F{x,y) ^ 



1 — X — y — xy 



whose coefficients fmn are called the Delannoy numbers, counts the number of lattice paths 
from (0,0) to {m,n) with allowable steps (1,0), (0, 1), and (1, 1). 
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To compute the asymptotics of the general diagonal coefficients fan,bn using the standard 
diagonal method, we fix x > small enough and try to find the small poles of 



1 r Fix^i^'/t^.t") 



dt 



27ri ./^_ t 

■dt 



1 r t'-^ 



2711 J^^ - x^l°- - t°-^^ - xV^i" 

Since ^ — x^l°- — t°-^^ — x^l°-t°- has a zero of multiplicity 6 at t = when x = 0, it follows 
that F{x^/'^ /t^ ^t"") /t has a single pole s of order h satisfying Y\tclx-,o s{x) = 0. Thus we 
have that 

G{x) = Res(F(x^/7^^^")/^;^ = s) 

= lim {{t - s)'F{x'/yt\ tn/t) , 

where is the derivative with respect to t. Patiently tracing through Leibniz's rule 
for the iterated derivative of a product, we could express this limit in terms of s and 
g{t) := (t — s)~''{t^ — x^l"" — t""^^ — and use it to find an algebraic equation satisfied 

by G. However, even with the help of Maple, this seems unlikely for general a and h. 

Thus, already at the first step of the standard diagonal method, it seems we are thwarted, 
the cause of the problem being the parameters a and h occurring in the exponent: the larger 
they get, the greater the difficulty in finding a defining equation for G. < 

So we should keep to main diagonal asymptotics when using the standard diagonal 
method. But even this poses problems in three or more variables, as illustrated by the 
next example. 

Example 2.3 {Ternary Words). The trivariate generating function 

1 



F{x,y,z) 



1 — X — y — z 

counts the the number of words fimn over a ternary alphabet, {0,1,2} say, that have / 
zeros, m ones, and n twos (an easy consequence of the symbolic method as described in 
|FSj ). An easy combinatorial argument shows that fimn = C^m'n"^)' which one could 
apply Stirling's formula and derive the asymptotics. 

To compute the asymptotics of the main diagonal coefficients fnnn using the standard 
diagonal method instead, we iterate the contour integration process. First 



1 f F{x,y/t,t)_^^ 



J J. t 

1 r 1 



■dt, 



27ri J^^ -t2 + (^i-x)t-y 

which has a single simple pole approaching zero as y and x approach 0, namely, s 
|(1 + X — Y^l — 2x + x^ — Ay). The residue is then 

G{x,y) := ^ . 

V 1 - 2x + x2 - 4y 
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Second, since 

lim G(x/t,t) = lim — , = , 

^-0 x^o ^i-2x/t + x^/t^ - it yr^' 

G{x/t,t) is a holomorphic function of t in some annulus about t = 0. Thus computing 

1 f G(x/M) ,, 1 / ' 



2iri J-,, t 2ni J.,^ t-/! - 2x/t + x'/f - 4t 

will give us the diagonal generating function. Proceeding, we employ Maple and find that 
two singularities of G{x/t,t)/t approach as x approaches 0, namely 



1 , , q(x) 1 1^/1,, q(x)\ 



24 

and its conjugate, where ^(x) = (— 36x+216x^ + l+24-v/— Sx^ + 81x^)^^^ and q{x) = x—^. 
While both of these singularities are algebraic in x, Maple can not compute the sum of 
their residues in reasonable time. 

The problem is that the residue sum, that is, the diagonal generating function ^ (n^n n)'''"' 
is not algebraic (see jSta99l Exercise solution 6.3] for instance). While it is D-finite — the 
diagonal of any rational function in any finite number of variables is D-finite |Lip88| — 
it not trivial to find a closed form or defining differential equation for it. (In the present 
case we obtain, using the Maple package Mgfun, the defining differential equation 27x(x — 
l)/"(x) + (54x — l)/'(x)+6/(x) = 0, which has a hyper geometric solution). Moreover, even 
given the defining differential equation, we still require the theory of univariate singularity 
analysis of D-finite functions to compute the asymptotics. But this theory has not been 
worked out in general. Indeed, as Philippe Flajolet has informed us, certain aspects of it, 
such as the so-called connection problem, might not be computable! <1 

3. A New Diagonal Method 

To transcend the limitations of the standard diagonal method, to go beyond main diag- 
onal asymptotics of bivariate generating functions, we take a multivariate approach. When 
the diagonal method works, it does at least produce an explicit algebraic formula for the 
diagonal generating function, which our method here does not. However, it is normally 
much more useful to have an asymptotic expression for the coefficients. In order to achieve 
this, we can apply a multivariate singularity analysis directly to F(x). Such an analysis 
has been recently developed by Baryshnikov, Pemantle, and Wilson in |P W02j . JPW04j . 
and [bp], the relevant parts of which we now summarize and adapt to our needso 

Let P C be the open domain of convergence of F and write -F(x) = /(x)/J(x) 
for some / and J holomorphic on an open domain D' containing the closure 2 ^) ^'^id 
relatively prime in the ring of holomorphic functions on D' . In all examples in this paper, 
the representation of -F as a quotient in fact holds on all of C^, but the extra generality is 
sometimes useful in applications. Let V be the complex variety {x G : J(x) = 0}. 

A critical point of F for n G (N^)'^ is a solution of 

J(x) = 

nrfXj Ji(x) = niXdJd(x) {i < d). 

Here Ui denotes component i of n, and Ji denotes the partial derivative of J with respect 
to component i of the domain of J, conventions we adopt for all vectors and functions 
throughout and which should cause no confusion in context. Let CRiT(n) denote the set 



-'^ As a whole, the multivariate singularity analysis developed in the aforementioned articles is much more 
general than what we present in this article, as it applies not only to rational functions but also to locally 
meromorphic functions. For a more complete summary including further examples, see the forthcoming 
article [PW] . 
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of all critical points of F for n. For generic directions n, this set is finite, being a zero- 
dimensional complex variety. The main situation in which CRiT(n) is infinite occurs when 
J defines a binomial variety {x | x** — x'^}, in which case CRlT(n) is empty for all but one 
direction and uncountable otherwise. Such examples can be analysed by a variant of the 
methods shown here. 

A contributing point of F for n is a critical point that influences the asymptotics of 
the coefficients of F in the direction of n. Let n denote the set of all such points. While n 
is ill-defined here, its functional role will become clear from the next two theorems. 

Theorem 3.1. If CRiT(n) is finite, then 

• CRiT(n) contains exactly one point, call it c, that lies in the positive orthant of 
M'^, and c e n; 

• all other members of n must lie on the same torus as c; 

• in the case where J = 1 — P for some aperiodic power series with nonnegative 
coefficients P, n = {c}. 

Here, the torus of a point c G is the set {x S : Vi < d \xi\ = |cj|}, and a power 
series is aperiodic if the Z-span of its monomial vectors is all of Z'^. 

A point c G V is a smooth point if V is a smooth complex manifold in a neighborhood 
of c, or equivalently, if Ji(c) / for some i (see |BK86l page 363] for instance). For 
simplicity of presentation, we deal only with smooth points in this article. This is the 
generic case, although interesting examples are not always generic. For more on the case 
of non-smooth points, see |PW04] . 

Theorem 3.2. Let n = (oin, . . . , a^n) for some oi, . . . , S N"^. If n consists of a single 
smooth point c such that CdJd{c) 7^ and q is a simple zero of Xd ^ J{ci, . . . , Cd-i,Xd), 
then 

\j=o 

for some bj G C as n ^ co. If n consists of a finite set of smooth points ci, . . . ,Cj. each 
satisfying the hypotheses above, then 

fn-- 

as n ^ 00. 

Theorem 3.3. In Theorem 13^ 




61 







/(c) 



-CdJd(.c)y/{2TTr-^h{J,c)' 

assuming h{J,c) ^ 0. Here h{J,c) is the determinant of the matrix H, whose entries are 

TJ^ ('^rnJl + Cd{JdJlm — JmJld — 'hJdm H Y^Jdd] 

'd^d \ 'Jd 



^d'^d 

o2 



H;; — — ^ -I- -2-72 ( J? + Cd{JdJll — '^Jl'Jld + -T-^dd)] 

CdJd cpi V Jd J 



'-d'^d 

where l,m < d and / 7^ m. 



Proof. The formula for bo comes from |PW021 Theorem 3.5]. (Please note the typo therein: 
the ZdHd should be a —ZdHd-) We prove the formula for H only. 
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For easy reading, let x = {xi, . . . Since J is holomorphic and Jrf(c) / 0, there 

exists a holomorphic function g guaranteed by the Implicit Function Theorem such that 
in some open ball around c, J(x, g(x)) = 0. As shown in |P W02j . H is the Hessian (matrix 
of second partial derivatives) of 

d-l 

m := +log g{cW'\. . .,Cd-ie''''-^) + i V - logg(c) 



evaluated at 9 = 0. 

Now, from the Implicit Function Theorem we also get 

J/(x,5r(x)) ,s 
D; = — -rj^ — 7^ (/ < d); 



Jdi'x,g{x)) 



J 11 



Dm Ji{^, = Jim + Jld Dm 5 = Jim " Jld-T {I < d,m< cl). 

Jd 

Thus, up to an additive constant, T>i f = —icie^^'g'^Ji/Jd, and for l,m < d with / / m we 
have 



=0 



-2t-. "^l , -iJd^m Jl - JlDm Jd 

-9 ^m9-j-+g ^2 

■Jd -Jd 



-iJl , 2f -2t-. '^l , -iJdDl Jl - JlBi Jd 
9=0 Jd V Jd Jd 



The result then follows by plugging in the Implicit Function Theorem equations above and 
simplifying. □ 

Several of our examples involve main diagonal asymptotics for symmetric J. In this case 
h{J,c) simplifies greatly. 

Proposition 3.4. If CRiT(n) is finite, n = (n, ...,n), J(x) is symmetric in x, and 
Jd{c) 7^ 0, where c is the contributing point that lies in the positive orthant of M'^, then 
c = (c, . . . , c) for some positive real c, and 

d-l 



h{J, c) = d 



l + -r{Jdd 
Jd 



Jdi) 



Proof. By Theorem 13.11 CRiT(n) contains exactly one point c that lies in the positive 
orthant of W^. By the symmetry of the critical equations J = and Xj Jj = XdJd < d) 
induced by the symmetry of n and J, any permutation of c's coordinates is also a critical 
point lying in the positive orthant of W^. Since there can be only one such point, c = 
(c, . . . , c) for some positive real c. 

Since J is symmetric and c = (c, . . . , c), all first partial derivatives of J are equal at c, 
all mixed second partial derivatives of J are equal at c, and all diagonal second partial 
derivatives of J are equal at c. Thus the entries of the matrix of Theorem 13.31 simplify to 



6=0 



e=o 



l + -r{Jdd 
■Jd 



J, 



dl, 



2Dimf 



8=0 



where l,m < d and I ^ m, and the determinant simplifies to h{J,c) = a"^ ^d, where 



l + i:iJdd 



Jdi) 



□ 



Our new diagonal method is simply the application of Theorems 13. H[3?2t and l3.3l directlv 
to F. We illustrate this now with several examples. 
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Example 3.5 {Zigzag-free Binary Words). Consider again the bivariate generating 
function 



^-^ 1 — X — y + xy — x^v^ 



m,n 



of Section [2l To compute the asymptotics of /„„ using the new diagonal method, let 
n = (n, n) and J = 1 — x — y + xy — x^y^. Then CRiT(n), the solution set of J = 
and nxJx = nyJy, comprises (1/0, 1/0), (-0, -0), (1 ± i\/3)/2, and (1 ± i\/3)/2, where 
= (1 + ^/5)/2. By Theorem EH c := (1/0,1/0) is a contributing point and the only 
such point, since none of the other critical points lie on the same torus as c. It is also a 
smooth point. By Proposition 13. 4t /i(J, c) = 4/(— 5 + 3\/5), and by Theorems 13.21 and \'6.'6\ 

f ^2n ^ 

Inn ~ , ■ 

This agrees with our answer in Section [2] since 

02 = 2/(3 - \/5). < 
Example 3.6 {Lattice Paths). Consider again the bivariate generating function 

Fix,y) ^ 



1 — X — y — xy 

of Section [2l To compute the asymptotics of fan,bn using the new diagonal method, let 
n = {an,bn) and J = l — x — y — xy. Then CRiT(n), the solution set of J = and 
bnxJx = anyJy, comprises (— where L = \/a'^ + b"^. By Theorem 13.11 c := 
(~^^' ~^a^) only contributing point. It is also a smooth point. By Theorem! 

HJ,c) = " . and by Theorem 



an,bn 



, r \ —an / T, I r \ ~bn / , 

-a + L\ I —b + L\ / ab 



L{a + b-L)^27m 



Here, where the standard diagonal method failed, the new diagonal method allowed us to 
compute the asymptotics quite easily (with the help of Maple). This success comes from 
a general attribute of the new diagonal method: the parameters a and b appear as factors 
instead of as powers in the equations we need to solve. < 

Example 3.7 {Ternary Words). Consider again the trivariate generating function 

1 



F{x,y,z) 



1 — x — y — z 

of Section [51 To compute the asymptotics of fan,bn,cn^ let n = {an, bn, cn) and J = 1 — 
X — y — z. Then CRiT(n), the solution set of J = 0, cnxJx = anzJz, and cnyjy = bnzJ^ 
comprises c := (a + 6+ c)^^(a, b, c). By Theorem 13. H c is the only contributing point. It is 
also a smooth point. By Proposition 13.41 /i(J, c) = ab{a + 6 + c)/c'^, and by Theorems 13.2 
and 13.31 

(a + 6 + c) ja + b + c 1 

J an,bn,cn ~ 



f^anjjbn^cn ]/ 27rn ' 

in agreement with Stirling's formula. 

Here, where the standard diagonal method failed because of the nonalgebraic character 
of the diagonal generating function, the new diagonal method gave the answer easily, even 
off the main diagonal. < 

Rather than multiply examples unnecessarily, we finish with an example involving an 
arbitrary number of variables. Many more examples of multivariate asymptotics using 
similar methods can be found in |PWj . 
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Example 3.8 (Alignments). In a problem from computational biology, one is given a 
finite set of strings of various lengths over some finite alphabet and needs to count the 
number of ways of inserting blanks into these strings to make them all of equal length. The 
strings represent genetic sequences, and inserting blanks represents aligning genetically 
similar segments among the sequences. 

More specifically and abstractly, a (d, ri, . . . , r^)- alignment is a table of d {0, l}-strings, 
one in each row, such that the sum of the entries of row j equals rj (1 < j < d) and 
no column contains all zeros. (The ones correspond to letters of the genetic sequence and 
the zeros to blanks.) As shown in ^GHOWQOj . the generating function for the number of 
(d, •)-alignments is 

To compute the asymptotics of fn,...,n (the case where all rows have the same sum; 
the genetic sequences all have the same length) using the new diagonal method, let n = 
(n, . . . ,n), / = 1/2 and J = 1 - (1/2) ULii'^ + ^i)- Then CRiT(n), the solution set of 
J = and nxiJi = nxd-Jd (* < d) comprises {2^/'^uji — 1, . . . ,2^/'^uJi — 1), where the oji 
{i < d) are the dth roots of unity. By Theorem 13. H c := (2^/'' — l,...,2i/'^-l) is the only 
contributing point. It is also a smooth point. By Proposition 13.41 h{ J,c) = d2^^~^^/'^, and 
by Theorems 13.21 and 13.31 

" ^ ^ ' (2i/<i - i)2{d2-i)/2d^d(^n)'^-i ' 

in agreement with [GHOW90] . (Please note the typo therein: the 2^^^~'^^/'^^ in the formula 
on page 139 should be a 2(i-'=')/2'=). 

A more important case for biological applications is to count alignments whose minimum 
block size is bounded below by a constant. A block of size 6 in a (d, •)-alignment is a d x 6 
submatrix of the alignment with contiguous columns, none of which contain a zero. As 
shown in [RT] . the generating function for (d, •)-alignments with blocks of size at least h is 



G{xi,...,Xd) 



l + {l-p){l-t + t^) -t'^ + t^+i 
I - t - [p - I - t){l - t + t^) 



where p = W^^ii^ + Xi) and t = ni=i ^i- Note that when 6 = 1, G = F, as expected. 

To compute the asymptotics of gn,...,n (which to our knowledge has only been done 
up till now for d = 2 via the standard diagonal method [RTj ) let n = (n, ...,n) and 
J = 1 — t — (p— 1 — t){l — t + t^). Then CRiT(n) comprises all vectors xl such that 
j{x) := J{xl) = (as seen from the equations nxjjj = nxdJd)- Let c := cl be the critical 
point in the positive orthant of M'' (so c is the unique positive solution of j{x) = 0) and 
note that c < 1, because j{x) < for x > 1. By Theorem 13.11 c S n and all other points 
of n lie on the same torus as c. Since 

1 + ^ 

G(xi, . . . - 



l-(p-l-t)(l + ^) 

on the open polydisk {x G : Vi \xi\ < 1}, since this polydisk contains the torus of c, 
and since (p — 1 — t)(l + Y~t) is an aperiodic power series with nonnegative coefficients, 
we can apply Theorem 13.11 to conclude that c is the only contributing point. It is also a 
smooth point, as is seen most easily from the form J = 1 — (p — 1 — t){l + Thus by 
Theorem 13.21 
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where 60 is the constant given in Theorem 13.31 Note that both c and 60 depend on b. < 
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